A^-body Efimov states from two-particle noise 
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The ground state energies of universal A^-body clusters tied to Efimov trimers, for A'^ even, are 
shown to be encapsulated in the statistical distribution of two particles interacting with a background 
auxiliary field at large Euclidean time when the interaction is tuned to the unitary point. Numerical 
evidence that this distribution is log-normal is presented, allowing one to predict the ground-state 
energies of the A''-body system. 



While noisy correlators are generally regarded as an 
impediment for Monte Carlo calculations, there is of- 
ten a physical mechanism underlying the appearance of 
noise. Understanding the form of the statistical distri- 
bution can lead not only to better methods for extract- 
ing reliable quantities from numerical calculations, but 
also to insight into the physics itself. In this work, we 
illustrate this principle using a correlator for two parti- 
cles with an interaction tuned to give a bound state at 
threshold, called the unitary point. From the distribution 
of this correlator we are able to extract the energies of 
2iV-body Efimov states, deeply bound universal systems 
of bosons or distinguishable fermions tuned to unitarity, 
which are tied to Efimov trimers [2 [5] ■ Efimov physics 
has enjoyed a resurgence of interest among the atomic, 
nuclear, and condensed matter communities due to ad- 
vances in ultracold atom experiments, particularly with 
recent experimental evidence for three- and four-body 
Efimov states displaying universal characteristics [3HH]j 
meaning their low-energy behavior is independent of the 
details of the interaction. However, theoretical informa- 
tion about the existence and properties of higher-body 
systems has been limited to that from direct measure- 
ments of TV-body states using non-perturbative numeri- 
cal methods [9Hl2]. 

Recent lattice studies of many-fermion systems at uni- 
tarity have shown that the correlators display distribu- 
tions with log- normal characteristics [131 E] . Using our 
two-body correlator we establish a deep connection be- 
tween Efimov physics and the log-normal distribution. 
We then present lattice data which indicates that the 
correlator is log-normal to within 2%, and use this knowl- 
edge to derive an analytical prediction for the energies of 
2iV-body Efimov states. Finally, we compare our results 
to those from numerical calculations. 

To begin, we will define the two-body correlation func- 
tion of interest and study its distribution by calculating 
the moments. The Lagrangian consists of two flavors of 
nonrelativistic particles [52] interacting through a point 
interaction with coupling n, 
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(1) 



Performing a Hubbard-Stratonovich transformation gives 
the Euclidean path integral in terms of an auxiliary field, 



K = {dr -y^/2M) + K(j). (2) 
The two-particle correlation function is given by 

C2(T)-(*(T)vl/t(0)), (3) 

where 



^'(t) = j dxidx2A{xi,X2)-4^^{xi,T)-4)\^{x2,T) (4) 

annihilates a two-body state, with wavefunction A at 
time T, which has non-zero overlap with the two-body 
ground state. Integrating out the ^ fields gives the cor- 
relation function as a path integral over (j) field configu- 
rations only. 



C2{T) = (C2(0,T))^ 



V(j){AeiKfS2{(l),T) 



2?0(detA')^e--'''^^'''^5'/>' , 



(5) 



where S2{4>,T) is the two-particle propagator from Eu- 
clidean time r = to T on a given background field, 0, 
and we have one power of det K for each flavor. Using 
open temporal boundary conditions, one may show that 
det is a constant independent of 4>, and therefore may 
be disregarded [15J. The use of open boundary condi- 
tions is justified so long as we restrict our arguments to 
zero temperature (large Euclidean time). 

By inserting a complete set of energy eigenstates, 
in Eq. [3| one may show that for large Euclidean time, T, 
C2 will be dominated by the ground state, 



C2(T)=^(M^(0)| 



T—foo 



,(6) 



(2) 

where is the ground state energy of the two-body 
system and Z2 gives the overlap of the operator "if with 
the ground state. 
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Following an argument of Lepage |16| . we determine 
the second moment of the correlator as follows, 

M2{T) ^ (|C2(^,T)|2)^ 

= ^ /l?(/.(52[<^,r])2e-/'^^'^'-5*'. (7) 
^0 J 

We see that this quantity is just the correlator for two 2- 
body propagators, corresponding to the four-particle sys- 
tem. Because there is no (anti-)symmetrization between 
the four (fermions)bosons, each particle must correspond 
to a different flavor. In order to properly describe a four- 
flavor theory we should also require four powers of det K 
in the partition function, Eq. [Sj However, because det K 
is a constant in our formulation, the partition function is 
unchanged when the number of flavors is increased. 
For the large Euclidean time limit we have, 

M2iT) Z,e-<'^^ (8) 

T— i-oo 

where Eq'^^ is the ground state energy of the four-particle, 
four-flavor system, and Z4 gives the overlap of ^I^^ with 
the four-particle ground state. We may generalize this 
argument for all moments, 

Mn{T) Z2Are-^o'"'^ , (9) 

T— >-oo 

where E^^'' is the ground state energy of the 2 A^-particle 
state in the 2iV-flavor theory and Z2N represents the 
overlap of '^^ with the 27V-particle ground state, which 
is assumed to be non-zero but may be arbitrarily small. 

Let us now consider the ground-state of these systems 
in the unitary limit, where the two-body ground state 
energy is zero by construction. The three-body system 
exhibits the well-known Efimov effect [U [5] , consisting of 
a series of bound trimers whose energies are separated by 
a factor of ^ 515. The ground-state is stabilized against 
collapse by an effective 3-body interaction which becomes 
relevant for low energy physics. Hence, the addition of 
a third particle to the conformally invariant two-particle 
system introduces an energy scale, A^;, to which the 3- 
body ground-state binding energy may be related, -Eq — 

Theoretical studies |9l HOl [T7H20] as well as recent ex- 
periments [3 IH] indicate that the spectrum of the four- 
body, four-flavor theory consists of a set of two bound 
tetramers tied to each Eflmov trimer. Provided the UV 
physics is fixed such that the ground state is sufficiently 
far from the cutoff, the ratio of the lowest four-body state 
to the lowest Efimov trimer has been shown to be a uni- 
versal constant [3 [101 [13 [HI HO] . This implies that there 
is no relevant four-body scale for this system, so A^; re- 
mains the only scale, giving Eq'^^ — —a^AE- 

One may postulate that no new relevant scales will 
emerge for the A^-body system at unitarity, and numeri- 
cal evidence for up to iV = 40 suggests that this is true 
[nilinidl]- Accordingly, the energy of the A^-body system 
will be Eq^^ = —gnAe- 



Returning to the probability distribution, we may com- 
bine the above relation between energies at unitarity with 
the moments of the distribution, giving 

7Wr'*(r)->Z2we°="'^^^ (10) 

Now we shall discuss the observed behavior of the dis- 
tribution using lattice results. In |13 it was shown that 
lattice calculations of strongly interacting non-relativistic 
systems tend to display heavy-tailed distributions at 
large Euclidean time. In particular, it was found that 
the distributions are approximately log-normal (LN), i.e. 
that the logarithm of the quantity of interest obeys the 
normal distribution. It was also shown that an expan- 
sion around LN, called the cumulant expansion, can be 
used to extract a reliable mean for the correlator. Using 
this method, the moments of In C are calculated and the 
correlator is given by InC(r) = limAr^_j.oo lnC^^"''(T), 
where 

lnC(^»)(r) = f]^^. (11) 

n=l ^' 

Here, k„(T) is the n-th cumulant of the distribution for 
In C (0, T) , and the expansion may be cut off for finite 
TVk provided the distribution is close enough to LN that 
the series converges. For the LN distribution, k„ = for 
n > 3, so the expansion may be cut off exactly at N^^ — 2. 

A histogram of the two-particle correlator at unitarity 
is shown in Fig. [l] for large Euclidean time. The distribu- 
tion of the logarithm of the correlator is also shown in the 
inset. The data was generated using the lattice method 
developed in [TS]. A set of momentum-dependent in- 
teractions are tuned to systematically correct for lattice 
artifacts, bringing us closer to the unitary point by ef- 
fectively setting the range of the interaction and the first 
few shape parameters to zero. Preliminary calculations 
have shown that an L = 16 box is sufficient to eliminate 
finite volume effects for the three-body state [21], which 
is the largest of the TV-body bound states and therefore 
most susceptible to finite volume errors. Finally, there 
is a hard momentum cutoff in our formulation. We find 
that in the absence of an explicit three-body interaction 
the resulting energy cutoff is approximately 335 times the 
three-particle ground-state energy per particle. 

We see that this correlator exhibits a very heavy tail, 
while the logarithm of the correlator appears to be Gaus- 
sian, characteristic of the LN distribution. To understand 
physically why this occurs, let us compare the moments 
derived above to those of the LN distribution. For the 
LN distribution the n-th moment is 

where fj,, a are the mean and standard deviation, respec- 
tively, of the logarithm of the correlator. 

We may extract values for these parameters by setting 
Eq. |9] equal to Eq. [T2] for n = TV = 1 and n = TV = 2 
to find n = \e''q^T , ct2 = -E^'^'r. Thus, for the special 
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FIG. 1: Histogram of the two-body correlator at large Eu- 
clidean time, T. Inset: Histogram of the logarithm of the 
correlator at the same value of T. 



(2) 

case of unitarity, where Eq — 0, we find that and 
a are no longer independent. This impHes that there is 
only one scale controlling all moments of the distribution. 
The moments become, 



1)E'„>T _ in(n-l)a4AET 



(13) 



A comparison of Eq. [TO] and Eq. [13] shows that the 
moments given by the LN distribution display the same 
scaling behavior with Ae,T as the moments predicted by 
the Efimov spectrum. Therefore, systems in which the 
Efimov effect plays a role are Ukely to have distributions 
with LN characteristics. In turn, having a LN distribu- 
tion implies an infinite tower of universal A^-body states 
whose energies are proportional to a single scale. Note 
that the LN distribution is the maximum entropy dis- 
tribution of its class: once the parameters fi and a have 
been fixed by the energies Eq^'' and Eq'^'' , the distribution 
should be LN in the absence of further constraints on the 
higher moments. 

Finally, assuming the distribution of the two-body cor- 
relator is exactly LN, we can deduce the energies of all 
2A^-body states by equating Eq. [9]with Eq. 13 

1 



E, 



{2N) 



-NiN^l)E', 



(4) 



(14) 



We see that the iV-dependence of the energy relation is 
equivalent to the number of pairwise interactions between 
dimers. 

For an approximately LN distribution one could imag- 
ine systematically improving this relation by numerically 
calculating third and higher cumulants of lnC2, using 
this information to correct the distribution, for example 
via the principle of maximum entropy, and recalculat- 
ing the moments of C2. We would like to emphasize 
that this approach does not simply translate a difficult 
many-body problem into an equally difficult problem of 
extracting large moments of a distribution, but rather 
small moments of In C2 ■ 

Remarkably, to within a few percent, we find that the 
distribution is indeed LN. In Fig. [2] (green points) we 



plot the cumulant expansion for the lattice data. Recall 
that for the LN distribution, this expansion converges at 
= 2. Thus, the discrepancy between the expansion 
cut off at Ni^ = 2 and the result after convergence may 
be used to quantify how close to LN a distribution is. 
We find that our lattice data is LN to within ^ 2%. This 
small discrepancy is likely due to the finite time extent 
used and sensitivity to lattice artifacts for large moments 
of C2. 

The energies for an exactly LN distribution (Eq. 14) 
are plotted in Fig. [3j along with results from a numeri- 
cal calculation employing a model potential. Comparing 
with [TD], we find agreement at the 10 — 30% level, with 
greater discrepancy for larger A^. 

The inset shows the results for Eq / E4 from three sepa- 
rate numerical calculations, [3], [in]j and [TT]. The three 
points from left to right represent a movement toward a 
more universal regime (see details in [101 E]). One sees 
a trend toward the LN result as the universal regime is 
reached. For larger A^, universality may be increasingly 
difficult to reach in numerical calculations as the more 
compact states begin to probe the details of the poten- 
tial model used in these calculations. 

One may ask by how much we can deform the relation 
between energies and still recover a distribution that is 
LN in appearance. In particular, could the small discrep- 
ancy from LN seen in the lattice data give the 10 — 30% 
shift in energies corresponding to those in [TOj ? 

To answer this question, we created mock distributions 
by expanding about a LN distribution and fitting the 
undetermined coefficients so that the moments gave the 
energies calculated in TO]. We created several of these 
distributions, using different parameterizations and fit- 
ting different combinations of energies. The results for 
the cumulant expansion, Eq. |11[ of one of the mock dis- 
tributions is plotted in Fig. [21 along with that for the 
lattice data of the two-body correlator. We find that the 
discrepancy for all mock distributions is ~ 17% — 30%, 
which is comparable to the difference in energies from the 
LN distribution and those of lOj. We conclude that the 
energies implied by our two-particle correlator should be 
those of Eq. [ill likely to within a few percent. 

The LN distribution implies other physical conse- 
quences that connect this distribution with universal be- 
havior. If we begin by choosing a wavefunction for our 
source that corresponds exactly to the solution for the 
two-body system, then Z2 = I and the constant of pro- 
portionality in Eq. 13 is, 



Z2N — -Z, 



iAr(A'-l) 



(15) 



We may interpret this relation as implying that the size 
of the system, and therefore the overlap, scales with A^ 
in the same way as the energy. 

Finally, note that if the two-body correlator is LN, 
then by extension the correlators for all 2Af-body Efimov 
states will also be LN, with rescaled parameters fi and 
a^. The logarithm of the sth moment for the 2A^-body 
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FIG. 2: (Co lor online) Deviation between the cumulant ex- 
pansion (Eq. 11 1 cut off at 7V„ and the converged value. The 



green points represent the lattice data for the two-body cor- 
relator tuned to unitarity on a 16^ x 1000 lattice, while the 
blue points are from a mock distribution created to reproduce 
moments corresponding to the energies found in [10|. 




FIG. 3; (Color online) Ratio of energies Ef'> /E'f' predicted 
by the LN distribution (green), compared to numerical calcu- 
lations, (blue) [10 , and (red) [TT]. Inset: Numerical results 
for £6/^4 from (black) g], (blue) [10], and (red) [11] • The 
solid line is the result predicted by the LN distribution. 



For odd numbers of particles, the correlator is not pos- 
itive so a LN distribution is not expected. However, in 
practice it is found that these distributions are approxi- 
mately LN with a small negative contribution. By fixing 
an additional ratio, such as E^/E^, one may extract ap- 
proximate relations between the energies for odd systems. 

To summarize, a connection between the log-normal 
distribution and Efimov physics has been established us- 
ing the distribution of the two-body correlator at unitar- 
ity. Using this connection, a novel method for obtaining 
the ground-state energies of 2A^-body Efimov states has 
been introduced. Because lattice data strongly indicates 
that the distribution of this correlator is LN to within a 
few percent, we calculate the expected ground-state en- 
ergies from the moments of the LN distribution. 

We note that the scaling of the energy per particle 
given by the LN distribution, Eq^^ /N ~ ^TV, implies 
that cutoff independence should not hold for arbitrarily 
large N . Thus, very large moments of the distribution are 
not expected to conform to those of the LN distribution. 
However, these moments correspond to the non-universal 
regime which is inherently of less interest than the uni- 
versal regime for which the LN distribution appears to 
be relevant. 

Given that the distribution is likely not exactly LN 
due to these large moments, analytical progress may be 
made by developing a perturbative expansion around LN, 
perhaps in the spirit of the semiclassical expansion intro- 
duced in [13]. Numerical efforts to reduce systematics in 
the lattice calculation of the moments of \nC2{(f>) would 
help to sort out true deviations from LN from lattice ar- 
tifacts, and any true deviations may be included in the 
overall distribution to obtain improved results for the en- 
ergies. 
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correlator is given by the energy relation, Eq. |14[ 



E, 



(2Ns) 



1 



T = ^Ns{Ns- l)ElpT . 



(16) 



Com par ing to the moments of the LN distribution 
(Eq. 



12|, we find = ^NE^o*h , = -N^E^^^'T 



r2 7p(4). 
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